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Abstract 

Many drugs exist in the crystalline solid state due to reasons of stability and ease of handling during the various stages of 
drug development. Crystalline solids can exist in the form of polymorphs, solvates or hydrates. Phase transitions such as 
polymorph intercon version, desolvation of solvate, formation of hydrate and conversion of crystalline to amorphous form 
may occur during various pharmaceutical processes, which may alter the dissolution rate and transport characteristics of the 
drug. Hence it is desirable to choose the most suitable and stable form of the drug in the initial stages of drug development. 
The current focus of research in the solid-state area is to understand the origins of polymorphism at the molecular level, and 
to predict and prepare the most stable polymorph of a drug. The recent advances in computational tools allow the prediction 
of possible polymorphs of the drug from its molecular structure. Sensitive analytical methods are being developed to 
understand the nature of polymorphism and to characterize the various crystalline forms of a drug in its dosage form. The 
aim of this review is to emphasize the recent advances made in the area of prediction and characterization of polymorphs and 
solvates, to address the current challenges faced by pharmaceutical scientists and to anticipate future developments. 
© 2001 Elsevier Science BY. All rights reserved. 

Keywords: Crystallinity; Polymorphs; Hydrates; Solvates; Formulation; Drug substance; Phase transformation; Characterization 



Contents 



1. Introduction 4 

2. Recent advances in the identification, prediction and characterization of polymorphs 6 

2.1. Types of polymorphism 6 

2.2. Packing polymorphism 7 

2.3. Conformational polymorphism 8 

2.4. Phase transformations in the solid state 9 

2.5. Prediction of polymorphs 11 

2.6. Directing the crystallization of specific polymorphs 12 

2.7. Characterization of polymorphs using a combination of analytical techniques 13 

3. Recent advances in the identification and characterization of hydrates and solvates 15 

3.1. Introduction to solvates and hydrates 15 

3.2. Structural aspects 15 



^Corresponding author. Tel.: + 1-612-6243-956; fax: + 1-612-6250-609. 
E-mail address: grant001@tc.umn.edu (DJ.W. Grant). 

0169-409X/01/$ - see front matter © 2001 Elsevier Science BY. All rights reserved. 
PII: S0169-409X(01)00097-7 




ELSEVIER Advanced Drug Delivery Reviews 48 (2001) 3-26 



4 



S.R. Vippagunta et al. I Advanced Drug Delivery Reviews 48 (2001) 3-26 



3.3. Phase transformation of hydrates and solvates 17 

3.4. Prediction of the formation of hydrates and solvates 18 

3.5. Characterization of hydrates and solvates 18 

4. Current challenges and future directions 18 

4. 1 . Origins of the challenges 18 

4.2. Phase transformations during processing 19 

4.3. Degree of cry stallinity 21 

4.4. Characterization of mixtures of polymorphs 21 

5. Conclusions 23 

References 24 



1. Introduction 

Most organic and inorganic compounds of phar- 
maceutical relevance can exist in one or more 
crystalline forms. When applied to solids, the adjec- 
tive, crystalline, implies an ideal crystal in which the 
structural units, termed unit cells, are repeated 
regularly and indefinitely in three dimensions in 
space. The unit cell has a definite orientation and 
shape defined by the translational vectors, a, b, and 
c, and hence has a definite volume, V, that contains 
the atoms and molecules necessary for generating the 
crystal. Each crystal can be classified as a member of 
one of seven possible crystal systems or crystal 
classes that are defined by the relationships between 
the individual dimensions, a, b, and c, of the unit cell 
and between the individual angles, a, f3, and y of the 
unit cell [1,2]. The structure of a given crystal may 
be assigned to one of the seven crystal systems, to 
one of the 14 Bravais lattices, and to one of the 230 
space groups [1]. All the 230 possible space groups, 
their symmetries, and the symmetries of their diffrac- 
tion patterns are compiled in the International Tables 
for Crystallography [3]. 

The common crystalline forms found for a given 
drug substance are polymorphs and solvates. Crys- 
talline polymorphs have the same chemical com- 
position but different internal crystal structures and, 
therefore, possess different physico-chemical prop- 
erties. The different crystal structures in polymorphs 
arise when the drug substance crystallizes in differ- 
ent crystal packing arrangements and /or different 
conformations. The occurrence of polymorphism is 
quite common among organic molecules, and a large 
number of polymorphic drug compounds have been 
noted and catalogued [4-7]. 

Solvates, also known as pseudopolymorphs, are 



crystalline solid adducts containing solvent mole- 
cules within the crystal structure, in either stoichio- 
metric or non stoichiometric proportions, giving rise 
to unique differences in the physical and pharma- 
ceutical properties of the drug. If the incorporated 
solvent is water, a solvate is termed a hydrate. 
Adducts frequently crystallize more easily because 
two molecules often can pack together with less 
difficulty than single molecules. While no definite 
explanations can be given, possible reasons include 
adduct symmetry, adduct-induced conformation 
changes, and the ability to form hydrogen bonds 
through the solvent molecules [2,8,9]. Desolvated 
solvates are produced when a solvate is desolvated 
and the crystal retains the structure of the solvate 
[10]. Desolvated solvates are less ordered than their 
crystalline counterparts and are difficult to character- 
ize, because analytical studies indicate that they are 
unsolvated materials (or anhydrous crystal forms) 
when, in fact, they have the structure of the solvated 
crystal form from which they were derived [11]. 

Because different crystalline polymorphs and sol- 
vates differ in crystal packing, and /or molecular 
conformation as well as in lattice energy and en- 
tropy, there are usually significant differences in 
their physical properties, such as density, hardness, 
tabletability, refractive index, melting point, enthalpy 
of fusion, vapor pressure, solubility, dissolution rate, 
other thermodynamic and kinetic properties and even 
color [12]. Differences in physical properties of 
various solid forms have an important effect on the 
processing of drug substances into drug products 
[13], while differences in solubility may have impli- 
cations on the absorption of the active drug from its 
dosage form [14], by affecting the dissolution rate 
and possibly the mass transport of the molecules. 
These concerns have led to an increased regulatory 
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interest in understanding the solid-state properties 
and behavior of drug substances. For approval of a 
new drug, the drug substance guideline of the US 
Food and Drug Administration (FDA) states that 
"appropriate" analytical procedures need to be used 
to detect polymorphs, hydrates and amorphous forms 
of the drug substance and also stresses the impor- 
tance of controlling the crystal form of the drug 
substance during the various stages of product 
development [11]. It is very important to control the 
crystal form of the drug during the various stages of 
drug development, because any phase change due to 
polymorph intercon versions, desolvation of solvates, 
formation of hydrates and change in the degree of 
crystallinity can alter the bioavailability of the drug. 
When going through a phase transition, a solid drug 
may undergo a change in its thermodynamic prop- 
erties, with consequent changes in its dissolution and 
transport characteristics [15]. 

Various pharmaceutical processes during drug 
development significantly influence the final crys- 
talline form of the drug in the dosage form. The 
various effects of pharmaceutical processing on drug 
polymorphs, solvates and phase transitions have been 
described in detail by Brittain and Fiese [16] and will 
be discussed in later chapters. Briefly, processes such 
as lyophilization and spray drying may lead to the 
formation of the amorphous form of drug, which 
tends to be less stable and more hygroscopic than the 
crystalline product. Also, processing stresses, such as 
drying, grinding, milling, wet granulation, oven 
drying and compaction, are reported to accelerate the 
phase transitions in pharmaceutical solids. The de- 
gree of polymorphic conversion will depend on the 
relative stability of the phases in question, and on the 
type and degree of mechanical processing applied. 
Keeping these factors in mind, it is desirable and 
usual to choose the most stable polymorphic form of 
the drug in the beginning and to control the crystal 
form and the distributions in size and shape of the 
drug crystals during the entire process of develop- 
ment. The presence of a metastable form during 
processing or in the final dosage form often leads to 
instability of drug release as a result of phase 
tran sf ormation [17]. 

Crystallization plays a critical role in controlling 
the crystalline form and the distribution in size and 
shape of the drug. The significance of crystallization 



mechanisms and kinetics in directing crystallization 
pathways of pharmaceutical solids and the factors 
affecting the formation of crystals have been re- 
viewed in detail by various researchers [12,18,19]. A 
crystalline phase is created as a consequence of 
molecular aggregation processes in solution that lead 
to the formation of nuclei, which achieve a certain 
size during the nucleation phase to enable growth 
into macroscopic crystals to take place during the 
growth phase. The factors affecting the rate and 
mechanisms by which crystals are formed are: 
solubility, supersaturation, rate at which supersatura- 
tion and desupersaturation occur, diffusivity, tem- 
perature, and the reactivity of surfaces towards 
nucleation. The various forces responsible for hold- 
ing the organic crystalline solids together, such as 
nonbonded interactions and hydrogen bonding, have 
been discussed in detail by Byrn et al. [2] and Etter 
[20]. 

Various analytical methods are being currently 
used to characterize the crystalline form of the drug 
during the various steps of processing and develop- 
ment. These methods have been reviewed recently in 
detail by many authors [7,10,21-25]. The single 
most valuable piece of information about the crys- 
talline solid, including the existence of polymorphs 
and solvates, is the molecular and crystalline struc- 
ture, which is determined by single-crystal X-ray 
diffractometry [2]. Powder X-ray diffractometry 
provides a "fingerprint" of the solid phase and may 
sometimes be used to determine crystal structure. 
Once the existence of polymorphism (or solvate 
formation) is definitely established by single-crystal 
and powder X-ray diffractometry, spectral methods, 
such as Fourier transform infrared absorption (FTIR) 
spectroscopy, Fourier transform Raman scattering 
(FT Raman) spectroscopy, solid-state nuclear mag- 
netic resonance (SSNMR) spectroscopy, ultraviolet 
and visible (UV-Vis) and /or fluorescence spectros- 
copy [23] may be employed for further characteriza- 
tion. Of special significance are thermal methods, 
such as differential scanning calorimetry (DSC), 
thermogravimetric analysis (TGA) and optical micro- 
scopy using a hot stage [24]. These methods are 
almost always employed for further characterization. 
Modulated (temperature) differential scanning 
calorimetry (MDSC) in combination with DSC and 
optical microscopy are able to identify the glass 
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transition of amorphous forms with much greater 
clarity and allow unique insights into the glass 
transitional and polymorphic behavior of drug sub- 
stances [26]. 

Because solid-state NMR spectroscopy can be 
used to study crystalline solids, as well as pharma- 
ceutical dosage forms, this powerful method is 
finding increasing application in deducing the nature 
of polymorphic variations [27], such as variations in 
hydrogen bonding network and molecular conforma- 
tions among polymorphs [28,29] and for the de- 
termination of molecular conformations and mobility 
of drugs in mixtures and dosage forms [2]. Solid- 

1 3 

state C-NMR in conjunction with the techniques, 
known as high power proton decoupling, cross 
polarization (CP), and magic-angle spinning (MAS) 
offers information not obtained readily by other 

1 3 

techniques. Recently, two-dimensional C- solid- 
state NMR spectroscopy has been used to study the 
three conformational polymorphs of 5-methyl-2-[(2- 
nitrophenyl)amino]-3-thiophenecarbonitrile [30]. Use 
of two-dimensional NMR and total suppression of 
spinning side bands (TOSS) pulse sequences allowed 
the separation of isotropic and anisotropic chemical 
shifts for the three forms. This is a very powerful 
method for analyzing differences in the chemical 
environment and is finding increased application in 
the study of conformational polymorphism. 

With advances in analytical methods, the current 
focus of research in the solid-state area is to under- 
stand polymorphism and pseudopolymorphism at the 
molecular level. Knowledge of the crystal packing 
arrangements and the various intermolecular forces 
involved in the different packing arrangements will 
help in the prediction and preparation of the most 
stable polymorphs of a given compound well in 
advance, to avoid surprises during product develop- 
ment. A current emphasis is on the development of 
software to predict crystal structures of polymorphs 
from molecular structures. A thorough understanding 
of the physicochemical properties of polymorphs and 
solvates (hydrates) is of primary importance to the 
selection of a suitable crystalline form and develop- 
ment of a successful pharmaceutical product. Bray et 
al. [31] have shown that, by thorough characteriza- 
tion of four different crystalline forms of L-738,167, 
a fibrinogen receptor antagonist by various analytical 



techniques, it was possible to determine the suitabili- 
ty of one or two forms for the development of 
pharmaceutical oral dosage forms. 

The present review aims to emphasize the recent 
advances made in the area of prediction and charac- 
terization of polymorphs and solvates, attempts to 
address the current challenges and problems faced by 
pharmaceutical scientists and intends to anticipate 
future development. This review does not attempt to 
provide solutions to the problems but attempts to 
review comprehensively the advances made in recent 
years to help address these problems. 



2. Recent advances in the identification, 
prediction and characterization of polymorphs 

2.1. Types of polymorphism 

Based on differences in the thermodynamic prop- 
erties, polymorphs are classified as either enantio- 
tropes or monotropes, depending upon whether one 
form can transform reversibly to another or not. In 
an enantiotropic system, a reversible transition be- 
tween polymorphs is possible at a definite transition 
temperature below the melting point. In a mono- 
tropic system, no reversible transition is observed 
between the polymorphs below the melting point. 
Four useful rules have been developed by Burger and 
Ramburger [32,33] to determine qualitatively the 
enantiotropic or monotropic nature of the relation- 
ship between polymorphs. These rules are the heat of 
transition rule, heat of fusion rule, infrared rule and 
density rule. 

If, by use of the above rules, it is established that 
the polymorphs of a particular drug are enantiotropic 
or monotropic, then the next goal is to define the 
thermodynamically stable (or metastable) domain of 
each crystalline phase of a substance as a function of 
temperature. The plot of the Gibbs free energy 
difference, AG, against the absolute temperature, T, 
gives the most complete and quantitative information 
on the stability relationship of polymorphs [22], with 
the most stable polymorph having the lowest Gibbs 
free energy. The AG between the polymorphs may 
be obtained using several techniques operating at 
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different temperatures, such as solubility [34] and 
intrinsic dissolution rate. Yu [35] has derived thermo- 
dynamic equations to calculate AG between two 
polymorphs and its temperature slope from the 
melting data. This method is essentially an extension 
of the heat of fusion rule, which is based on 
statistical mechanics. Extrapolating AG to zero gives 
an estimate of the transition temperature, from which 
the existence of monotropy or enantiotropy is in- 
ferred. The integration of different types of data 
provides the AG vs. T curve over a wide temperature 
range and allows the consistency between techniques 
to be checked [22] . Another approach to establish the 
order of stability among various polymorphs has 
been studied using pressure versus temperature plots, 
e.g., for sulfanilamide and piracetam [36]. This 
approach is based upon Ostwald's principle of least 
vapor pressure, according to which the stable poly- 
morph exhibits the lowest vapor pressure. The 
accuracy of this approach to establish the stability 
hierarchy among the polymorphs has been shown to 
be very much dependent on the accuracy of the 
experimental data. 

In recent years, the main focus of research has 
been the characterization of polymorphs arising from 
structural differences in the crystal lattice. It has 
been established for some time that organic mole- 
cules are capable of forming different crystal lattices 
through two different mechanisms. One of the mech- 
anisms is termed packing polymorphism, and repre- 
sents instances where conformationally relatively 
rigid molecules can be assembled into different 
three-dimensional structures through the invocation 
of different intermolecular mechanisms. The other 
mechanism is termed conformational polymorphism 
and arises when a nonconformationally rigid mole- 
cule can be folded into different arrangements, which 
subsequently can be packed into alternative crystal 
structures. The distinction between packing poly- 
morphism and conformational polymorphism is 
somewhat artificial because different packing ar- 
rangements impose different conformations on the 
molecules, however slight, and different conforma- 
tions will inevitably pack differently. The structural 
aspects associated with polymorphs have been re- 
viewed recently [2], as have the analogous features 
of solvate and hydrate systems [9]. In the next 



section, the results of some more recent investiga- 
tions are discussed. 

2.2. Packing polymorphism 

An investigation into the structures and charge 
densities of two polymorphs of p-nitrophenol has 
been performed with the aim of deducing the differ- 
ent modes of inter-molecular hydrogen bonding that 
lead to the formation of the two structures shown in 
Fig. la and b [37]. A detailed analysis of the charge 
density of the two forms indicates charge migration 
from the benzene ring region to the nitro and 
hydroxyl groups that accompanies the transformation 
of one form into the other. In addition, polarization 
of the oxygen lone-pair electrons was found to be 
substantially larger in the crystal forms than in the 
free molecule, resulting in considerably larger dipole 
moments in the solid state. 

During the study of a new crystal form (form I) of 




Fig. 1. Molecular packing diagrams of the (a) (3 polymorph of 
/7-nitrophenol, (b) a polymorph of /?-nitrophenol, showing 50% 
probability displacement ellipsoids ([37], reproduced with the 
permission of the American Chemical Society). 
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chlordiazepoxide, it was found that the heat of 
transition between the two forms (forms I and II) is 
rather modest, and kinetic factors permit the exist- 
ence of the metastable phase [38]. Both structures 
contain four crystallographically independent mole- 
cules linked in dimers through hydrogen bonding, 
but the dimers are packed differently to yield the two 
crystal forms. Because the dimers in the fundamental 
units are spaced differently in the two forms, it was 
proposed that the solid-state enantiotropic trans- 
formation entailed rearrangement of the dimer units. 

A different approach has been taken during an 
evaluation of the different structures formed by 
sulfathiazole [39]. Using a graph set approach to 
classify the known structural differences and simi- 
larities among the various forms, it became possible 
to identify packing motifs common to three of the 
four crystal structures. Fig. 2 shows the unit cells of 
the polymorphs I, II, III and IV, where molecules are 
paired as hydrogen-bonded dimers. At the end of the 
process, the authors were able to deduce possible 
links between the observed patterns of hydrogen 
bonding, processes of nucleation, and the crystal 
growth observed from a number of solvent systems. 
Interestingly, the analysis did not indicate a relation- 
ship between the appearance of a particular poly- 
morph from solution and the growth of its fastest 




Fig. 2. Unit cells of four polymorphs (I, II, III and IV) of 
sulfathiazole showing hydrogen bonds, with the dimer structure 
clearly discernible ([39], reproduced with the permission of the 
Royal Society of Chemistry). 



growing surface. Rather, it appeared as if the differ- 
ent solvents affected the process of polymorph 
formation through their effects on nucleation of the 
various forms. 

2.3. Conformational polymorphism 

The conformational polymorphism of the two 
forms of piroxicam pivalate has been studied in 
detail [40]. This compound is distinctive in that the 
high-melting form (polymorph 1) contains an un- 
anticipated array of associated molecules bound as 
centrosymmetric dimers through hydrogen bonding, 
with the amido nitrogen atom acting as the donor and 
the pyridine nitrogen as the acceptor (Scheme 1, 
structure I). The low-melting form (polymorph 2) 
contains molecules of two distinct conformational 
states coexisting in the same crystal (Fig. 3), but 
linked through different hydrogen bonding arrange- 
ments. This latter finding represents another unusual 
aspect of the crystallography of the substance. 

The inclusion of different solvent molecules in a 
crystal lattice can lead to the existence of different 
packing patterns, and has also been found to in- 
fluence the molecular conformation of paroxetine 
hydrochloride in two solvate forms [41]. One form 




ill 

Scheme 1. Molecular structure of piroxicam pivalate (I) [40], 
5-methyl-2- [2-(nitrophenyl)amino]-3-thiophenecarbonitrile (II) 
[30], 2'-acetamido-4'-[A^,A^-bis(2-methylcarbonylethyl)amino]-4- 
nitroazobenzene (III) [48]. 
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Fig. 3. Conformations of the two independent molecules of 
piroxicam pivalate (I) in polymorph 2. Thermal ellipsoids are 
drawn at the 40% probability level, and H atoms are shown as 
spheres of arbitrary size ([40], reproduced with the permission of 
the American Pharmaceutical Association). 



hemihydrate phase. Crystals of the isopropanol sol- 
vate decomposes in the open air at room tempera- 
ture, because the isopropanol molecules are released 
easily through the channel. The hemihydrate is 
relatively stable. 

In an impressive fundamental study, the poly- 
morphism of 5-methyl-2-[2-(nitrophenyl)amino]-3- 
thiophenecarbonitrile (Scheme 1, structure II) has 
been catalogued [42,43] and discussed in detail [2]. 
This compound was crystallized as six solvent-free 
polymorphs, each of which differed in the mode of 
packing and in molecular conformation. The differ- 
ent conformers yielded sufficient perturbations on the 
respective molecular orbital so that a variety of 
crystal colors (red, orange, and yellow) were ob- 
served. To obtain a more detailed evaluation of the 
relative stability, the authors considered a partition- 
ing of polymorphic energy differences into lattice 
and conformational contributions, and were able to 
deduce general trends that appeared valid in the 
absence of hydrogen bonding. The act of crystalliza- 
tion was found to feature an interplay of opposing 
forces, with perpendicular molecular conformations 
being favored in fluid solutions, while a preference 
for planar /high dipole conformers existed in most 
crystal forms, as shown in Fig. 4 [42]. Hie unusual 
polymorphism displayed by this system may result 
from one or more of the following factors: the 
preference for perpendicular conformations in solu- 
tions, the preference for planar /high dipole confor- 
mers in crystals, the formation of inter- and in- 
tramolecular hydrogen bonds, and the thermody- 
namic tendency towards low energy and high en- 
tropy. 



was obtained as a hemihydrate, and the other as the 
solvate of isopropanol (2-propanol). In the unit cell 
of the hemihydrate, one finds two protonated parox- 
etine and two chloride ions together with one water 
molecule. Interestingly, the two paroxetine molecules 
are conformationally nonequivalent, and exhibit a 
number of different bond angles and torsion angles. 
In the other form, the unit cell contains one proton- 
ated paroxetine molecule, one chloride ion, and one 
isopropanol molecule disordered along a molecular 
channel. Furthermore, the conformation of the parox- 
etine molecule in the isopropanol solvate is different 
from either molecular conformation observed in the 



2.4. Phase transformations in the solid state 

Studies of phase transformations in the solid state 
are important, because the sudden appearance or 
disappearance of a crystalline form can threaten 
process development, and can lead to serious phar- 
maceutical consequences if the transformation occurs 
in the dosage forms. Hence, an understanding of the 
kinetics and mechanism of phase transformations is 
of practical importance. The rearrangement of mole- 
cules into a new structure during phase transforma- 
tion may or may not involve a solvent or vapor 
phase. To explain the mechanism of solid- solid 
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Fig. 4. Comparison of conformational energies and crystal ener- 
gies of the various polymorphs of 5-methyl-2-[2-(nitro- 
phenyl)amino]-3-thiophenecarbonitrile (II). Form R (red prisms, 
m.p. 106.2°C), form Y (yellow prisms, m.p. 109.8°C), form OP 
(orange plates, m.p. 112.7°C), form ON (orange needles, m.p. 
114.8°C), form YN (yellow needles, m.p. not measurable) ([42], 
reproduced with the permission of the American Chemical Socie- 
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physical transition, four steps have been proposed: 
(a) molecular loosening in the initial phase; (b) 
formation of an intermediate solid solution; (c) 
nucleation of the new solid phase and (d) growth of 
the new phase [2]. In an interesting study, 
Skwierczynski [44] has proposed a two-environment 
model to describe the decomposition reaction kinet- 
ics of a crystalline solid, aspartame. The decomposi- 
tion reaction of aspartame is a simple unimolecular 
thermally-induced aminolysis and the reaction 
proceeds under anhydrous conditions, i.e., water is 
not a reactant [45] . This model links the chemistry of 
the solid-state reaction with the molecular mobility 
of the reactant as the reaction proceeds. The advan- 
tage of this model is that it can be used to determine 
the shelf life of a product from kinetic data gathered 
at elevated temperatures. Apart from solid- solid 
physical transformations, solution-mediated physical 
transformations among polymorphs are also known 
to occur in processes, such as wet granulation and 
during dissolution testing. 

While the majority of studies have probed the 
equilibrium properties of polymorphic and solid-state 
solvated systems, relatively few have been concerned 
with the dynamics of phase transformation. Byrn et 



al. [2] have reviewed briefly the aspect of poly- 
morphic interconversions and the factors affecting 
the transformations. In one study, the contribution of 
hydrogen bonding to the a — > (3 phase change of 
resorcinol has been detailed [46]. The a form is 
more stable than the p form at room temperature, but 
is less dense than the (3 form. The transition of 
a — » |3 at an estimated transition temperature of 
337 ±1 K is accompanied by an increase in crystal 
density, with the structure shifting from an open 
array of molecules (linked through hydrogen bond- 
ing) to a denser structure resembling molecular 
crystals. Through the use of a simple potential 
model, it was concluded that, during the phase 
transformation, the energy of the hydrogen bonds 
decreases along with the extent of such bonding. The 
energy liberated by this process is almost offset by 
the enhanced Van der Waals energies associated with 
the increase in crystal density, and consequently the 
transition enthalpy is rather small. Accompanying 
the shifts in hydrogen bonding is a number of 
effective proton transfers, altering the covalent and 
ionic portions of the crystal. It was also learned that 
the increase in entropy produced from the redistribu- 
tion of protons was of the same order of magnitude 
as the entropy of the phase transition. 

A number of spectroscopic techniques have been 
used to study the processes associated with a poly- 
morphic transition of 2-(2,4-dinitrobenzyl)-3- 
methylpyridine [47]. The two interconverting struc- 
tures coexisted over a temperature range of at least 
8-9°C. The phase change was associated with a 
molecular tautomerization that translated through the 
collective changes of a large number of molecules, 
yielding domains having definite short-range order. 
The slowly evolving spectroscopy that took place 
above the transition temperature was interpreted as 
the annealing of domains into a long-range ordered 
system. The process of phase transformation ap- 
peared to consist of an initial fast redistribution of 
the mole ratio of the coexisting phases, followed by 
a much slower process involving a macroscopic 
relaxation of the system. Although local thermo- 
dynamic equilibrium was thought to exist in in- 
dividual domains, the magnitude noted for the tem- 
perature range of the phase transition was proposed 
to arise from nonequilibrium conditions existing 
among the various types of domain. 
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A combination of solid-state 15 N-NMR spectros- 
copy and X-ray crystallography was used to study 
polymorphic transitions in an azobenzene dyestuff, 
2' -acetamido-4' - [7V,iV-bis(2-methylcar- 
bonylethyl)amino]-4-nitroazobenzene (Scheme 1, 
structure III) [48]. This work established that the 
structure of one polymorph was disordered, and that 
the process of phase transformation entailed a crank- 
shaft-type motion of the azo linkage. The ORTEP 
plots of the two molecular conformations for the 
X-ray structure determination of structure III at 293 
K are shown in Fig. 5. Selective polarization inver- 
sion and band shape-fitting experiments were used to 
deduce the thermodynamic parameters of the ex- 
change process. 

Raman spectroscopy was used to study the effect 
of pressure on the phase transitions in hexa- 
methylbenzene and hexa(methyl-d 3 )benzene [49]. 
The form 11^ form III transition of the partially 
deuterated substance was found to take place at a 
lower pressure relative to that of the analogous 
hexamethylbenzene compound, which was attributed 
to differences in the energies of the intramolecular 
methyl torsional vibration in the two crystal forms. 
In another study performed by the same group, the 
effects of both temperature and pressure on the phase 
transitions of tetrafluoro- 1 ,4-benzoquinone were con- 
sidered [50]. In this system, the changes in en- 




Fig. 5. ORTEP plots of the two molecular conformations (1 and 
2) for the X-ray structure determination of 2 '-acetamido-4 '- [N,N- 
bis(2-methoxycarbonylethyl)amino]-4-nitroazobenzene (Scheme I 
Structure) at 293 K. Thermal ellipsoids are shown at 30% for 
clarity, with conformer 2 being represented by the solid lines and 
Conformer I by the dotted lines ([48], reproduced with the 
permission of the American Chemical Society). 



vironmental conditions were found to influence a 
number of intermolecular and intramolecular vibra- 
tional modes, yielding conformational changes that 
in turn produced the observed phase transitions. 

2.5. Prediction of polymorphs 

The main challenge in managing the phenomenon 
of multiple solid forms of a drug is the inability to 
predict the number of forms that can be expected in a 
given case. This prediction would involve quantifica- 
tion of the myriad intermolecular forces within any 
proposed crystal structure as well as the ability to 
postulate the likely packing modes for a given 
molecule in all its configurations [10]. Accurate 
theoretical prediction of polymorphs from studies of 
molecular dynamics and crystal structure generation 
would be of outstanding importance in drug research 
[36]. 

More research is now being directed towards 
developing computational tools to understand the 
nature of polymorphism and to predict polymorphic 
forms at an early stage in the drug development 
process. The recent developments in computational 
chemistry allow the prediction of possible poly- 
morphic forms based only on the molecular structure 
of the drug. The Polymorph Predictor, from Molecu- 
lar Simulations, is currently the only commercial 
software package that can predict the possible poly- 
morphs of an organic compound from its molecular 
structure [51]. The package developed by Karfunkel 
and co-workers [52-54] uses a Monte Carlo simu- 
lated annealing approach to generate thousands of 
possible crystal packing alternatives for a given 
molecule. Each of the unique crystal structures is 
then subjected to a lattice energy minimization to 
obtain the relative stability ranking of the various 
packing possibilities and the resulting lowest-energy 
structures are the potential polymorphs. This method 
has been successfully employed to generate known 
polymorphs of primidone (Fig. 6 A and B) and 
progesterone, starting from the molecular structures 
alone [55]. It has also been used to predict poly- 
morphs for a range of small molecules and to predict 
unknown polymorphic structures of 4-amidinoin- 
danone guanylhydrazone, a selective inhibitor of S- 
adenosylmethionine decarboxylase [56], and of as- 
pirin [57]. 
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Fig. 6. Comparison of the crystal structure of primidone A (A) 
versus the most likely packing arrangement, frame 7 (B) ([55], 
reproduced with the permission of Elsevier Science). 



The theoretical predictions of lattice energies, 
entropies, morphologies and polymorphs should 
stimulate experimental activities and vice versa. The 
current crystal-modeling efforts have the potential of 
producing more quantitative tools for bridging struc- 
tures and properties, which could help in creating 
solid forms with desired properties [22]. There are 
many limitations in using computational methods for 
predicting polymorphs theoretically. The first limita- 
tion is that the ab initio screening is useful only for 
nonionic rigid molecules. For more complex sys- 
tems, the method is very useful for generating 
plausible crystal structures, but it is not accurate 
enough to determine which of these possible struc- 
tures can actually be crystallized [58]. In addition, 
the limitations in computer power can restrict the use 
of this method for predicting polymorphs of complex 
molecules. An issue of concern is that the existing 
methods only predict the lattice energies, which 
relate to internal energies or enthalpies of the 
crystals. However, the relative thermodynamic 
stability of polymorphs is determined by the Gibbs 
free energy, which is a linear function of both 
enthalpy and entropy. Predictions of the relative 
stability of polymorphs will be more accurate when 
the entropies, as well as lattice energies, are consid- 
ered. Application of molecular dynamics may enable 
the entropies to be calculated. Hence, no general 
method is currently available for the prediction or 
interpretation of the properties of complicated poly- 
morphic or pseudopolymorphic systems. 



2.6. Directing the crystallization of specific 
polymorphs 

Complementing the different computational meth- 
ods for predicting the stable polymorphs of a given 
compound, various experimental methods are also 
being employed extensively to control the type of 
polymorph formed during the crystallization process. 
Many studies have reported the role of additives in 
controlling the outcome of the crystallization pro- 
cess. Some of the preselected additives are capable 
of inhibiting the nucleation and /or growth of the 
unwanted polymorphs. For the first time, the role of 
reaction by-products in controlling polymorph ap- 
pearance of a drug has been reported [59] . This drug 
is sulfathiazole that is known to exist as polymorphs, 
forms I, II, III and IV, that differ in the hydrogen- 
bond network. Form I was found to be different from 
the other three forms as a result of a different 
hydrogen bonding at the aniline moiety of the 
molecule. From studies of the hydrogen-bonding 
pattern, it was predicted that the ethamido derivative 
of sulfathiazole could selectively control the forma- 
tion of form I over other forms by entering the 
growing face of form I without disrupting the 
structure (Fig. 7a). Because a similar effect was not 
possible with the other forms, incorporation of the 
ethamido derivative in the other forms should inhibit 
their growth (Fig. 7b). Experimentally, it was shown 
that the ethamido by-product stabilized form I over 
the other polymorphs. This study clearly shows that 
the combination of crystal morphology and the 
hydrogen-bond network analysis of the different 
polymorphs offer a new and powerful approach to 
understanding and controlling polymorph appearance 
and stability in the presence of additives. 

A similar approach was also applied to stabilize a 
metastable a conformational polymorph of l- 
glutamic acid using additives [60]. Methods such as 
DREIDING and TRIPOS force fields were used to 
select appropriate additives which could mimic the a 
and |3 conformations. Four additives were chosen for 
this study of which two were present exclusively in 
the p conformation and theoretically should selec- 
tively inhibit the crystallization of the p phase and 
thus stabilize the metastable a phase. Experimental- 
ly, it was proven that the additives, by virtue of their 
conformation, were able to selectively inhibit the 
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Fig. 7. Possible binding interaction of ethamidosulfamide in the fastest growing faces of (a) form I and (b) form IV of sulfathiazole ([59], 
reproduced with the permission of Elsevier Science). 



appearance of the stable p polymorph of L-glutamic 
acid by interfering with either the nucleation rates or 
the growth rates and thus stabilize the metastable 
form. These studies demonstrate clearly that the 
molecular packing and intermolecular hydrogen 
bonds are the main features, which make possible the 
conformational discrimination. The use of conforma- 
tional mimicry to stabilize the metastable structures 
of conformational polymorphs now offers a powerful 
tool for the prediction and development of robust 
processes for the control of polymorphic systems. 

2. 7. Characterization of polymorphs using a 
combination of analytical techniques 

The common techniques often fail to differentiate 
definitively between two structurally similar poly- 
morphs. Hence more advanced techniques or a 
combination of techniques need to be used to avoid 
errors of interpretation and in the identification of 
polymorphs [24]. Combinations of techniques are 
being employed currently for the characterization of 
crystalline pharmaceutical solids. For example, con- 
ventional single-crystal X-ray diffractometry and 
polarized microscopy were of no use in distinguish- 
ing the two forms I and II of roxifiban, a very 
promising cardiovascular drug, because of the rela- 
tively small crystallite sizes of the polymorphs. 
Hence, transmission electron microscopy (TEM) and 



synchrotron X-ray diffraction techniques were em- 
ployed to characterize the unit cells of the two forms. 
By coupling the highly resolved synchrotron powder 
X-ray diffraction data shown in Fig. 8, with in- 




Q.A. 1 

Fig. 8. A partial comparison of (a) a synchrotron pattern of 
polymorph II of roxifiban, collected using a wavelength of 

o 

1.00006 A with (b) a conventional X-ray diffraction pattern using 
CuKct radiation in a region where there are many overlapping 
peaks. The patterns are plotted as a function of Q = 2w/d = 4tt 
sin Qlk to remove the effects of different wavelengths ([61], 
reproduced with the permission of the American Pharmaceutical 
Association). 
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formation obtained from TEM diffraction patterns, 
the unit cell parameters of the two forms of roxifiban 
were determined [61]. Similarly, the three modi- 
fications, I, II and III, of the nonsteroidal antiinflam- 
matory drug, tiaprofenic acid, could not be dis- 
tinguished by the two traditional spectroscopic meth- 
ods, FTIR and FT-Raman spectroscopy. The modi- 
fications can only be distinguished by a combination 
of thermoanalytical and powder X-ray diffractomet- 
ric methods [62]. 

Another example, to which a combination of 
techniques has been successfully applied to identify 
the various conformational polymorphs of a drug, is 
the characterization of the solid forms of neotame 
[29]. Neotame, A^-(3,3-dimethylbutyl)-L-aspartyl-L- 
phenylalanine methyl ester, a new high-potency 
sweetener exists in the following phase-pure crys- 
talline forms: monohydrate, the most stable crys- 
talline form of neotame under ambient conditions, a 
methanol + water solvate [63], a methanol solvate 
[64], an amorphous anhydrate [29] and a crystalline 
anhydrate (form A; [65]). The authors conducted a 
systematic study of the conversion of the monohy- 
drate under vacuum to a mixture of anhydrate forms 
followed by the reconversion of the anhydrate to the 
monohydrate upon exposure to moisture under am- 
bient conditions. No significant changes were ob- 
served in the powder X-ray diffraction patterns 
during part of the reconversion process, suggesting 
that no change in lattice structure had occurred. 
However, the solid-state 13 C-CP-MAS NMR spectra, 
indicated the presence of several forms of neotame 
during the reconversion (Fig. 9). This discrepancy in 
the results between the two techniques was attributed 
to the conformational change of neotame molecules 
during reconversion, without significant change in 
unit cell parameters. This example indicates that both 
solid-state 13 C-CP-MAS NMR spectroscopy and 
powder X-ray diffractometry are needed to analyze 
mixtures of solid forms of conformationally flexible 
molecules, such as neotame. 

1 3 

A combination of solid-state C-NMR spectros- 
copy and single crystal X-ray diffractometry also has 
been used to examine the solid-state tautomerism of 
acetohexamide [66,67]. Polymorphism of the anti- 
diabetic drug acetohexamide has been investigated 
by numerous techniques. On the basis of FTIR data, 
form A of acetohexamide has been proposed to exist 




ppm 

Fig. 9. Resonance signal of the phenyl carbon (C-15) attached to 
the side chain in the 13 C-CP-MAS NMR spectra of neotame 
anhydrate: (a) sample generated by placing the original monohy- 
drate under vacuum (~1 Torr) for 3 days; (b-e) sample after 
being sealed in a jar for 2, 4, 6 and 8 days, respectively; (f) 
sample after being exposed to a relative humidity environment of 
84% for 12 days ([29], reproduced with the permission of the 
American Chemical Society). 



in the enol-tautomeric state, whereas form B has 
been proposed to be in the keto-tautomeric state. 
Using NMR and crystal structure data it was firmly 
established that both these acetohexamide polymor- 
phic forms are present in the keto-form. Hence the 
combination of solid-state NMR spectroscopy and 
X-ray crystallography provided strong evidence that 
both forms of acetohexamide exist in the keto- 
tautomeric state and are truly polymorphic. 



S.R. Vippagunta et al. I Advanced Drug Delivery Reviews 48 (2001) 3-26 



15 



3. Recent advances in the identification and 
characterization of hydrates and solvates 

3.1. Introduction to solvates and hydrates 

It has been estimated that approximately one-third 
of the pharmaceutically active substances are capable 
of forming crystalline hydrates [68]. The water 
molecule, because of its small size, can easily fill 
structural voids and because of its multidirectional 
hydrogen bonding capability, is also ideal for linking 
a majority of drug molecules into stable crystal 
structures [2]. The mere presence of water in a 
system is not a sufficient reason to expect hydrate 
formation, because some compounds, though they 
are soluble in water, do not form hydrates. It is the 
activity of water in the medium that determines 
whether a given hydrate structure will form. Solvates 
may be formed when a pure organic solvent or a 
mixture of solvents is used as the solvent for 
crystallizing the compound. Guillory [69] has dis- 
cussed the various methods of preparation of hy- 
drates and solvates in detail. Because solvates be- 
have similarly to hydrates, common analytical tech- 
niques can be used for characterization of solvates 
and hydrates. 

3.2. Structural aspects 

Crystalline hydrates, based on their structure may 
be classified into three categories. The first category 
(class I) are the isolated site hydrates, where the 
water molecules are isolated from direct contact with 
other water molecules by intervening drug mole- 
cules, e.g., cephradine dihydrate. The second cate- 
gory (class 2) are channel hydrates where the water 
molecules included in the lattice lie next to other 
water molecules of adjoining unit cells along an axis 
of the lattice, forming channels through the crystal, 
e.g., ampicillin trihydrate. Hie channel hydrates can 
be subclassified into two subcategories. One category 
comprises the expanded-channel or nonstoichiomet- 
ric hydrates, which may take up additional moisture 
in the channels when exposed to high humidity and 
for which the crystal lattice may expand or contract 
as the hydration or dehydration proceeds effecting 
changes in the dimensions of unit cells, e.g. 
cromolyn sodium. The other subcategory comprises 



the planar hydrates, which are channel hydrates in 
which water is localized in a two-dimensional order, 
or plane, e.g., sodium ibuprofen. The third category 
(class 3) of crystalline hydrates are the ion-associated 
hydrates, in which the metal ions are coordinated 
with water, e.g., calteridol calcium [8,9]. 

In this section, some examples of nonstoichiomet- 
ric hydrates and their characterization will be dis- 
cussed in detail because these forms pose a special 
challenge in dosage form development due to unpre- 
dictability of water content in the crystals. Following 
the work of Cox et al. [70] the unusual water uptake 
and formation of nonstoichiometric hydrates of 
cromolyn sodium was reinvestigated using single 
crystal X-ray diffractometry, PXRD, as well as by 
molecular modeling [71]. Cromolyn sodium, an 
antiasthmatic drug, exists as two liquid crystalline 
phases and a crystalline hydrate phase that sorbs and 
liberates water continuously and reversibly to give a 
continuous range of nonstoichiometric hydrates [70]. 
The changes in the PXRD patterns of the crystalline 
hydrate phase of cromolyn sodium in response to the 
surrounding relative humidity (RH) were explained 
in the light of the molecular and crystal structure of 
cromolyn sodium. Single crystal X-ray diffrac- 
tometry indicated the space group for cromolyn 
sodium as PI, a chiral space group, even though the 
molecule itself is achiral. The crystal structure of 
cromolyn sodium with five or six water molecules 
per cromolyn sodium molecule, solved at room 
temperature by Hamodrakas et al. [72], revealed the 
positions of only one sodium ion and two water 
molecules and showed that the second sodium ion 
and the other water molecules are disordered. Re- 
cently, the single crystal structure of cromolyn 
sodium at 76% RH, with 6.44 molecules of water 
was solved at 173 K by Chen et al. [71]. This work 
showed that the second undetermined sodium ion is 
disordered over three sites and that four of the eight 
water positions are partially occupied. Comparison 
of the crystal structures determined by Hamodrakas 
et al. [72] and Chen et al. [71] indicated that the 
cromolyn anion is flexible. In particular, the bond 
and torsional angles of the 2-hydroxypropane linking 
the two cyclic moieties, changed to accommodate 
lattice expansion or contraction resulting from water 
sorption and desorption by the crystals. As water is 
taken up, the relative occupancies of the sites of the 
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second sodium ion and that of water molecules 
change. As a result, the triclinic structure with a > 
90° approaches the monoclinic form with a ~ 90°. 
To summarize, the presence of large water channels, 
the flexibility of the 2-hydroxypropane link, the 
disorder of the second sodium ion (Fig. 10) and the 
disorder of the surrounding water molecules in the 
crystal lattice explain the reversible and nonstoich- 
iometric water sorption and desorption by cromolyn 
sodium. This study emphasizes the importance of the 
detailed single crystal structure in explaining many 
unusual physico-chemical properties of drug hy- 
drates. 

The muscarinic agonist, LY297802 tartarate {i.e., 
( + )-3-[3-(butylthio)-l,2,5-thiadiazol-4-yl]-l- 
azabicyclo[2.2.2]octane monohydrogentartrate}, was 
also found to exhibit an unusual tendency to form 
nonstoichiometric hydrates of variable, but specific 
composition, ranging from 0 to 0.5 mol of water 

1 3 

[73]. Solid-state C-NMR spectroscopy, in conjunc- 
tion with moisture sorption analysis and X-ray 
crystallography was used to provide unique insights 
into nonstoichiometric moisture sorption behavior. 
The PXRD patterns of the drug exposed to different 
RH values (0 to 75%), indicated neither a peak shift 
nor the presence of any new peaks, suggesting that 



the anhydrous and the hydrated forms of the drug are 
isomorphic. Fig. 1 1 shows the significant changes in 
the SSNMR peaks on exposure to different relative 
humidities and temperature, indicating that water 
incorporated into the crystal lattice changes the local 
chemical environment and causes the observed NMR 
changes. The incorporation of water into the crystal 
lattice of the drug was also confirmed by X-ray 
crystallography. The considerable hygroscopic ity of 
the drug was rationalized in terms of the similar 
crystal structures of the hydrated and nonhydrated 
forms, and hence no significant structural modi- 
fications are needed for the reabsorption of water 
into the solids. The rates of dehydration and rehydra- 
tion are largely determined by the size of the water 
channels and the strength of the hydrogen-bonding 
interactions that bind the water molecules in the 
channels. 

Another interesting study with different solvated 
forms of L-lysine monohydrochloride (LH) was 
conducted by Bandyopadhyay et al. [74]. LH was 
found to form: a pure methanol solvate at water 
activity, a w < 0.34, with methanol activity, a m > 0.7; 
a dihydrate at a w > ~ 0.65 with a m <0.45; and 
mixed solvates at intermediate values of a w and a m . 
It was found that the dihydrate and the mono- 
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Fig. 10. In hydrated cromolyn sodium, coordination environment of: (a) the first (ordered) sodium ion, Na(l), shown in two neighboring 
unit cells (A and B); and the second (disordered) sodium ion at the three partially occupied sites, (b) Na(2), (c) Na(3), and (d) Na(4). The 
striped circles represent the sodium sites. The open circles represent the oxygen atoms coordinated to Na(l). The dotted (gray) circles 
represent the oxygen atoms coordinated to Na(2), Na(3), or Na(4) ([71], reproduced with the permission of the American Pharmaceutical 
Association). 
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Fig. 11. Solid-state 13 C-NMR spectra of LY297802 tartrate {( + )- 
3-[3-(butylthio)-l,2,5-thiadiazol-4-yl]-l-azabicyclo[2.2.2]octane 
monohydrogentartrate} after storage at 0% humidity and before 
and after drying at elevated temperatures ([73], reproduced with 
the permission of the American Pharmaceutical Association). 

methanol solvated forms of LH possessed similar 
crystal structures, similar PXRD patterns, but dif- 
fered in the crystal habit. The crystal structures of 
LH hydrate and methanol solvate indicate that one 
molecule of methanol in the methanol solvate oc- 
cupies approximately the same volume as the two 
water molecules in the dihydrate. During dehydration 
in the presence of methanol in the crystallization 
medium, the loss of one or more water molecules 
from the crystal lattice was compensated by the 
gradual uptake of methanol into the crystal structure 
to satisfy the hydrogen bonding pattern within the 
lattice with minor rearrangements, giving rise to 
mixed solvates. The similarities of the crystal lattices 
of the dihydrate and monomethanol solvate explain 
the similarity of the PXRD patterns. 



3.3. Phase transformation of hydrates and solvates 

Phase changes due to hydration /dehydration and 
solvation /desolvation of pharmaceutical compounds 
during processing or in the final product may result 
in an unstable system that would effect the bioavail- 
ability of drug from solid dosage forms. Various 
types of phase changes are possible in solid-state 
hydrated or solvated systems in response to changes 
in environmental conditions, such as relative humidi- 
ty, temperature and pressure. For example, some 
hydrated compounds may convert to an amorphous 
phase upon dehydration and some may convert from 
a lower to a higher state of hydration yielding forms 
with lower solubility. Alternatively, a kinetically 
favored but thermodynamically unstable form may 
be converted during pharmaceutical processing to a 
more stable and less soluble form [8]. The phase 
transitions in hydrates and solvates can occur at 
various stages of dosage form development. Morris 
[9] has discussed the behavior of hydrates during 
processing, handling and storage of formulations in 
detail. 

The phase transformations associated with expo- 
sure to water, such as during solubility measure- 
ments, wet granulation processes, dissolution studies 
and accelerated stability tests are likely to occur via 
solution mediation. Solution mediated phase trans- 
formations depend upon the solution phase to pro- 
vide the mobility necessary to rearrange in the most 
stable form and hence are much faster than solid- 
state transformations. The rate of a solution-mediated 
transformation is proportional to the solubility of the 
species involved. Temperature, pressure and relative 
humidity may increase the rate of phase transforma- 
tion of hydrates by inducing mobility in the system. 

Solution-mediated phase transformations have 
been reported for many hydrate systems, such as 
theophylline crystals [17], eprosartan mesylate [75] 
and nedocromil sodium [76]. Ghosh and Grant [77] 
have addressed a common problem associated with 
the characterization of solvates which centers around 
the determination of solubilities of solvates and of 
nonsolvates that undergo phase transformation in the 
presence of an interacting solvent, such as solvation 
of nonsolvates in the solvent of crystallization or the 
desolvation of solvates in water. A thermodynamic 
cycle analogous to Hess's law but based on free 
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energies has been developed to predict the theoretical 
solubilities of l,2-dialkyl-3-hydroxy-4-pyridones, 
which form 1:1 formic acid solvates in the presence 
of formic acid, and of the 1:1 formic acid solvates 
which produce the corresponding unsolvated com- 
pounds in the presence of water. A good correlation 
was obtained between the solubility values measured 
by the standard extrapolation method and that calcu- 
lated by means of the thermodynamic cycle. 

Apart from identifying and characterizing the 
phases during various stages of drug development, it 
is very important to gain an understanding of the 
dehydration /hydration mechanisms and kinetics. 
Many models have been developed to account for the 
dehydration kinetics of the crystalline hydrates [78]. 
Nucleation is the most significant phenomenon in 
determining the transformation kinetics, that is, the 
rate of formation of a new phase [8]. The dehydra- 
tion kinetics to some extent will also depend upon 
the class of the hydrate system to which the drug 
belongs, particle size and morphology. Hie practical 
applications of understanding the dehydration kinet- 
ics, as indicated by Morris [9], are mainly the 
determination of the conditions for allowable expo- 
sure of bulk drug substances during development and 
processing, proper packaging, allowable temperature 
ranges for shipping, storage, and labeling of the final 
product, and the initial selection of a form for 
development. 

3.4. Prediction of the formation of hydrates and 
solvates 

Predicting the formation of solvates or hydrates of 
a compound and the number of molecules of water 
or solvent incorporated into the crystal lattice of a 
compound is complex and difficult. Each solid 
compound responds uniquely to the possible forma- 
tion of solvates or hydrates and hence generalizations 
cannot be made for a series of related compounds. 
Certain molecular shapes and features favor the 
formation of crystals without solvent; these com- 
pounds tend to be stabilized by efficient packing of 
molecules in the crystal lattice, whereas other crystal 
forms are more stable in the presence of water 
and /or solvents. There may be too many possibilities 
so that no computer programs are currently available 



for predicting the crystal structures of hydrates and 
solvates. 

3.5. Characterization of hydrates and solvates 

The common methods for the characterization of 
hydrates and solvates are polarized light microscopy 
and hot stage microscopy, DSC, TGA, Karl Fischer 
titrimetry, single-crystal X-ray diffractometry, pow- 
der X-ray diffractometry, and infrared spectroscopy. 
These methods have been reviewed in detail [21] and 
will also be discussed in detail in later chapters. 

Pressure DSC is gaining increasing popularity in 
the study of solvates and hydrates where dehydration 
reactions occur above or near the boiling point of 
water. Using conventional DSC, it is very difficult to 
measure the heats of dehydration and heat of vapor- 
ization separately, but if one conducts DSC experi- 
ments at elevated pressures, the two processes may 
be completely separated. The advantage of using 
pressure DSC is that the pressure can be precisely 
controlled and the solids can be subjected to a 
controlled temperature program while under substan- 
tially elevated temperatures. The influence of ele- 
vated pressures on the solid-state behavior of car- 
bamazepine dihydrate was studied by Han and 
Suryanarayanan [79]. In Fig. 12 it is shown that 
pressure DSC can separate the dehydration and 
vaporization endotherms of carbamazepine dihydrate 
during its conversion to the anhydrate form. Also the 
technique permitted the water liberated on dehydra- 
tion to remain in intimate contact with the anhydrous 
phase formed which could significantly influence its 
solid-state properties. 

The combined physical analytical techniques of 
thermogravimetry and infrared spectroscopy (TG/ 
IR) can permit identification of the solvent incorpo- 
rated into the crystal lattice. This combined tech- 
nique has been used to study formulated products, 
such as capsules and tablets [80]. 

4. Current challenges and future directions 

4.1. Origins of the challenges 

A series of flow charts and decision trees have 
been presented and discussed [11,22] that can be 
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Fig. 12. Differential scanning calorimetry (DSC) curves of car- 
bamazepine dihydrate at different pressures: (a) at atmospheric 
pressure in a conventional DSC cell, (b) at atmospheric pressure 
in a pressure DSC (PDSC) cell, (c) 100 p.s.i., (d) 200 p.s.i., (e) 
300 p.s.i., (f) 400 p.s.i., and (g) 600 p.s.i. (1 p.s.i. = 6.9 kPa) 
([79], reproduced with the permission of Elsevier Science). 



used by investigators to characterize the polymorphs 
and solvates of compounds under development or for 
registration with regulatory authorities. Due to the 
complex and nonconventional behavior of various 
organic drug molecules, there are many opportunities 
for research and development in the area of charac- 
terization of polymorphs and solvates. Some of the 
problems which are commonly encountered during 
characterization of crystalline solids and which need 
to be addressed are: disorder in the crystal lattice due 
to pharmaceutical processing leading to conversion 
of a crystalline phase to an amorphous material or 
phase conversion from one form to the other; 
quantitating the amount of single polymorph in a 
mixture of polymorphs; identifying the solid form of 
the active ingredient in the formulated product, 
particularly when the drug is a minor component in 
the presence of numerous other materials (excipi- 



ents); and the issues of disappearing polymorphs and 
the appearance of new polymorphs. In the following 
sections we will address some of these issues and 
some of the studies that have addressed these 
problems. 

4.2. Phase transformations during processing 

The effects of pharmaceutical processing on the 
crystalline state of drug polymorphs and solvates 
have been discussed recently by Brittain and Fiese 
[16]. Exposure to changes in temperature, pressure, 
relative humidity and comminution are encountered 
during processes such as drying, granulation, milling 
and compression. The stresses applied to crystals 
during pharmaceutical processing can cause defects 
in their crystal lattices, and contribute to lattice 
disorder, thus affecting the physical properties of the 
resulting powder [81]. This problem has been dis- 
cussed in detail by Byrn et al. [10]. Arising from 
different degrees of crystalline disorder, the difficulty 
in reproducing materials with the same properties is 
a major concern in the pharmaceutical industry. 

Milling, the last processing step in the production 
of bulk drug substance to reduce particle size, is 
often accompanied by a decrease in crystallinity due 
to the creation of lattice defects, beginning at the 
surface. The defects created by mechanical activation 
of the solid on the surface can migrate, transform, 
and change their number and nature. If the defects in 
the mechanically activated crystal heal to produce a 
crystal lattice different from the initial lattice, then a 
polymorphic transformation has taken place. Milling- 
induced polymorphic changes have been observed 
for many small drug molecules, such as fostedil, 
chloramphenicol palmitate, indomethacin and 
phenylbutazone [16]. Polymorphic transformation of 
the dipeptide sweetener, aspartame hemihydrate, can 
occur during milling [82]. Polymorph II of aspartame 
hemihydrate was found to transform to form I during 
ball-milling or on heating for 30 min at 160°C in the 
presence of steam as shown in the X-ray diffraction 
pattern (Fig. 13). The susceptibility of form II of 
aspartame hemihydrate to transform to form I has 
been attributed to the less symmetric crystal structure 
of form II compared to that of form I as studied by 
spectroscopic methods. 

Some authors, such as Huttenrauch [81] have 
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Fig. 13. Powder X-ray diffraction patterns of aspartame 
hemihydrate: (a) theoretical powder pattern calculated from the 
crystal structure of aspartame hemihydrate (form I) determined by 
Hatada et al., (1985) [99], (b) experimental powder pattern of the 
ball- milled aspartame hemihydrate (form I), (c) experimental 
powder pattern of aspartame hemihydrate that had been heated for 
30 min at 160°C in the presence of steam (form I), (d) experimen- 
tal powder pattern of aspartame hemihydrate after compression at 
250 MPa for 1 min (form II), and (e) aspartame hemihydrate as 
received (form II) ([82], reproduced with the permission of the 
American Pharmaceutical Association). 

suggested that the trauma to crystals during grinding 
may lead to a decrease in crystallinity, which should 
improve the compression capacity and dissolution 
rate of the drug molecules. This hypothesis was 
tested by studying the morphology, crystalline state, 
compression capacity and dissolution properties of 
native and ground crystals of aspirin and lactose 
monohydrate [83]. No significant increase in com- 
pression capacity was observed when native and 
ground crystals were compared. Only a slight in- 
crease in the dissolution rate was observed for 
ground aspirin crystals which was attributed to 
surface defects due to grinding resulting in improved 
crystal wetting. 

The effect of roller compaction on lattice defects 
and phase change has been examined for aspirin 
[84]. The water vapor sorption isotherm obtained for 
aspirin on roller compaction indicated much more 



water uptake than had been reported previously for 
other crystalline samples of aspirin under similar 
conditions. Various possibilities were suggested for 
the unusual water uptake of aspirin on roller compac- 
tion, such as formation of a polymorphic form of 
aspirin with a much greater affinity for water, 
formation of a crystalline hydrate in the aspirin 
sample, or significant reduction in particle size of the 
aspirin particles thereby providing an increased 
specific surface area for water vapor adsorption. It is 
also possible that roller compaction disrupted the 
crystalline order of some part of the aspirin crystals 
forming amorphous regions, which then take up 
relatively large quantities of water into their bulk 
structure. This example clearly indicates that pro- 
cesses, such as roller compaction, can introduce 
considerable disorder in the surface of crystals 
leading to a marked increase in the tendency to sorb 
water vapor. 

Thermal activation, like mechanical activation 
during processing, also results in a high-energy state 
of crystals that may reorganize into a different lattice 
arrangement resulting in a phase change. The ther- 
mal stability of drug substances is important, because 
formulations are often dried at elevated temperatures 
after wet granulations so that the tablets may contain 
small regions of high temperature (hot spots) during 
compression. Various examples have shown that a 
change of temperature may influence the stability of 
drug molecules [16]. The effect of low temperatures, 
such as during freeze-drying, on the crystalline form 
of the drug has also been studied. The formation of a 
new mannitol hydrate during freeze-drying has been 
reported [85]. The formation of a crystalline hydrate 
by an excipient during freeze-drying may have 
several practical consequences, such that the difficul- 
ty of removing bound water from the crystal lattice 
can significantly limit the drying rate, while the 
residual water that is not removed by freeze-drying 
may be a potential threat to product stability if it is 
released during storage. The mannitol hydrate 
formed during freeze-drying survived the typical 
drying cycle and converted to the anhydrous poly- 
morph of mannitol upon heating. 

Spray drying has also been shown to lead to loss 
of crystallinity in materials, by a combination of 
processes involving rapid solidification of dissolved 
material and solid-state transitions due to milling 
effects in the atomiser. Spray drying leads to conver- 
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sion of a crystalline phase to an amorphous state and, 
because the amorphous state is metastable with 
respect to the crystalline form, phase transformations 
are likely to occur within the shelf life of the 
pharmaceutical product, resulting in loss of quality 
and potency in the product [86]. 

In view of the significant effects that the state of 
disorder in crystalline solids caused by pharmaceu- 
tical processing can have on the properties of 
pharmaceutical solids, it is important to be able to 
assess the extent of disorder in a solid quantitatively 
down to very low levels. Various methods have been 
used to measure the percent disorder, such as using 
predetermined mixtures, measurements of X-ray 
powder diffraction, density and heats of crystalliza- 
tion which revealed limits of detectability down to 
about 10%. Using water vapor sorption measure- 
ments under very carefully controlled conditions, it 
was possible to detect disorder as low as 1% in 
milled samples of sucrose [87]. A comparison of the 
four methods mentioned above for estimating the 
percent disorder of milled samples of sucrose gave 
very consistent results, once the underlying factors 
that make these techniques sensitive to the con- 
centration of amorphous content were recognized 
and taken into account. 

4.3. Degree of crystallinity 

The previous section has emphasized that many 
pharmaceutical processes lead to a decrease in 
crystallinity of drug phases. Various studies have 
concluded that the formation of amorphous material 
during processing is highly undesirable. The amor- 
phous material, being in a thermodynamically meta- 
stable state, is susceptible to reconversion to the 
crystalline state, affecting many physico-chemical 
characteristics of the drug. A later chapter provides 
detailed coverage of amorphous materials. An esti- 
mation of the degree of crystallinity of a sample 
before and after processing poses one of the larger 
challenges facing the pharmaceutical field. Powder 
X-ray diffractometry is still the commonly used 
method for determining the degree of crystallinity, 
though this method suffers from some limitations 
due to peak broadening, amorphous halo, and pre- 
ferred orientation which make interpretation and 
quantitation difficult. DSC may not be a sensitive 
method for measuring crystallinity due to crystalliza- 



tion of the amorphous content at elevated tempera- 
tures and the effects of differences in heat capacity. 
Solution calorimetry has been proposed as an accur- 
ate method for analysis of percent crystallinity 
[11,88,89]. A decrease in the endothermic enthalpy 
of solution indicates a decrease in the crystallinity of 
the sample. However, differences in surface area 
produced by grinding or by other processing tech- 
niques can also result in changes in the heat of 
wetting of a sample. Judicious choice of solvent can 
be employed to reduce such surface effects, which 
themselves contribute to the observed crystallinity of 
the sample. 

Near infrared (NIR) spectroscopy is another tech- 
nique being used to measure the degree of crys- 
tallinity, and has also proved useful in studies of the 
polymorphism and water content of sugars. The NIR 
spectrum of a sample contains both physical and 
chemical information. Being noninvasive, nondes- 
tructive and operable at room temperature, the 
method is a valuable tool with which to assess 
changes in the amorphous and crystalline state of 
lactose [90]. NIR has been used to follow the 
changes in the amorphous state, the onset of crys- 
tallization, and the changes between a- and |3 -lac- 
tose, which accompany the onset of crystallization. 
In another study, the nucleation and crystallization 
kinetics of amorphous lactose was investigated by 
gravimetry in an automated vacuum moisture bal- 
ance. The combination of isothermal and nonisother- 
mal activation energies allowed the investigation of 
both crystal growth and nucleation mechanisms and 
led to the separation of activation energies for 
nucleation and growth [91]. 

4.4. Characterization of mixtures of polymorphs 

Another common problem encountered during 
drug development is quantitative control of the 
proportion of polymorphic forms present in a mix- 
ture. According to the US FDA regulations, the 
method of analysis for the proportion of forms must 
be validated, and also the proportion of forms must 
remain within stated limits through the retest date of 
the drug substance and potentially throughout the 
shelf life of the product. This is a very onerous 
requirement, especially if the forms have a tendency 
to interconvert. Byrn et al. [11] suggested that the 
best way to deal with this problem is probably by 
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developing methods to prepare only one crystal form 
and maintaining this form throughout processing. 
Powder X-ray diffractometry is often a useful meth- 
od to determine the percentages of polymorphs in a 
mixture. However, the detection limit is variable 
from case to case, and is sometimes as high as 15%. 
It is therefore important to develop sensitive ana- 
lytical methods with a lower limit of detection. 

Attenuated total reflectance (ATR) FTIR spec- 
troscopy has been shown to be valuable for the 
quantitative analysis of the polymorphic content of 
bulk pharmaceutical materials. The feasibility of 
using ATR-FTIR for the qualitative and quantitative 
analysis of mixtures of pharmaceutical polymorphs 
has been studied using three known polymorphs of 
ganciclovir as a model compound [92]. Definitive 
identification and quantitation of all three poly- 
morphs could be achieved using ATR-FTIR spec- 
troscopy in conjunction with partial least-squares 
modeling. This technique has many advantages, such 
as speed, nondestructiveness, relative ease of use, 
and most important, no sample pretreatment before 
measurement. 

Raman spectroscopy is another technique that is 
being widely used to quantitatively estimate the 
percentage of one polymorph in a mixture of poly- 
morphs. FT-Raman spectrometry offers many advan- 
tages, the most prominent being, minimal sample 
preparation, sensitivity to polymorphism, and nonin- 
terference from water. Two polymorphs of 
fluconazole were characterized using FT-Raman 
spectroscopy and principal components regression 
using cross validation provided quantitative analysis 
of the percentage of one polymorphic form in the 
mixture of other forms [93]. A novel sample holder 
was developed whereby the sample is held in an 
NMR tube which is rotated around its axis and at the 
same time moved up and down. This method of 
sample presentation leads to a large increase in the 
volume studied and is important for inhomogeneous 
samples for which sub-sampling is a problem. Pos- 
sible degradation of the sample through heating by 
the laser can also be avoided [94]. 

X-Ray powder diffractometry is still the common 
method for the quantitative estimation of polymorphs 
in a mixture of polymorphs. This method requires 
that at least one high-intensity peak unique to each 
form is available for intensity measurements and that 



the plot of the peak intensity ratio as a function of 
the weight ratio of the components should result in a 
straight line. Modern computer controlled X-ray 
powder diffractometers now permit quantitative anal- 
ysis of multicomponent mixtures using the complete 
powder diffraction profile rather than a limited 
amount of low-angle integrated intensity data. Artifi- 
cial neural networks (ANNs) in quantitative X-ray 
powder diffractometry were used successfully to 
identify and quantify the two known modifications of 
ranitidine hydrochloride even when the weight frac- 
tion of one polymorph in the mixture was as low as 
0.01 [95]. ANNs have been used mostly in problems 
of pattern recognition and modeling, and is therefore 
useful in deciphering the pattern in diffraction data 
from polymorphic mixtures. The ANNs model pre- 
dicted concentration precisely, accurately, and with 
minimal bias through a wide range of ratios of the 
two known ranitidine hydrochloride polymorphic 
forms in a mixture (Fig. 14). This method minimizes 
the problems associated with preferred orientation 
and overlapping X-ray lines. The same group of 
researchers has shown the potential of ANNs in 
combination with DRIFT spectroscopy to analyze the 
polymorphic purity of crystalline ranitidine hydro- 
chloride as a bulk drug and as an active ingredient of 
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Fig. 14. Predicted concentrations of the two polymorphic forms 
(form 1 and 2) of ranitidine hydrochloride by the response surface 
methodology (RSM), a statistical modeling method, and by the 
artificial neural networks (ANNs) method, plotted against mea- 
sured percentage of form 2 ([95], reproduced with the permission 
of Elsevier Science). 
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a tablet formulation. Simultaneous identification and 
quantification of all the ingredients in the tablet 
formulation was possible. This study has shown that 
the complex problem of quantifying a drug in 
mixtures containing two or more components with 
overlapping spectra can be solved by the DRIFT- 
ANNs technique [96]. 

Another technique, which is gaining popularity in 
the quantitative analysis of mixtures of polymorphs, 
is solid-state 13 C-CP-MAS NMR spectroscopy. This 
method was preferred for quantitative analysis of 
polymorphic mixtures of the herbicide, pen- 
dimethalin, which exists as two polymorphs with 

1 3 

different colors and crystal habits [27]. C-NMR 
provided the most sensitive and definitive evidence 
of the transition from the yellow to the orange form. 
This method enabled as little as 2% of orange 
pendimethalin to be determined in a sample consist- 
ing mostly of the yellow polymorph. It is the least 
invasive of the instrumental methods and can be used 
to detect the ratio of the two polymorphs in solid 
formulations. 

A related challenge faced by the pharmaceutical 
industry is the determination of the polymorphic 
nature of the drug in the presence of excipients in a 
dosage form especially when the active drug is 
present as a low percentage of the overall formula- 
tion [97]. This problem can be addressed by develop- 
ing sensitive techniques with lower limits of de- 
tection or by using a combination of techniques. 

5. Conclusions 

In order to save time and cost it is very important 
to choose the most suitable form of the crystalline 
drug in the initial stages of drug development. In 
recent years a good deal of research has been 
directed towards achieving this goal. Systematic 
isolation and early characterization of the largest 
number of possible forms of a drug reduces the 
chances of surprises at the late production stage due 
to identification of a new crystalline form or phase 
change. With the development of more sophisticated 
computational tools, the main focus of many inves- 
tigators is to be able to predict all the possible forms 
of a drug from its molecular structure. Understanding 
the origins of the multiple solid forms of a drug 



molecule, either due to differences in packing ar- 
rangement or conformation of the molecules, be- 
comes the first step in prediction. Single crystal 
X-ray diffractometry and solid-state NMR spectros- 
copy are two techniques that are gaining increased 
application in determining the various crystal struc- 
tures and the origins of polymorphism and pseudo- 
polymorphism of a particular drug. When crystal 
structures can be calculated with certainty, it will be 
possible to predict the various polymorphs of a 
compound and this information could be used to 
guide experimental studies. This goal may be dif- 
ficult to achieve owing to the complex molecular 
structures of new organic molecules and the presence 
of several molecules in each asymmetric unit, but the 
future development of improved force fields and 
increased computational speeds, may make it achiev- 
able. 

Improved experimental methods leading to more 
accurate and detailed phase diagrams are also finding 
increased use in determining the stability of various 
polymorphs. It is important to make every effort to 
prepare and to identify the most stable polymorph in 
order to guide the selection of the optimal form for 
development. The emergence of sensitive methods 
and the use of combination techniques, facilitate the 
identification and the more accurate characterization 
of the various polymorphs of a drug molecule. One 
of the main analytical challenges faced by the 
pharmaceutical analyst is the development of better 
quantitative methods for identifying a single poly- 
morph in a mixture of polymorphs and for determin- 
ing the percentages of amorphous or crystalline 
content of the drug. More and more sensitive meth- 
ods are being developed to address this problem. 

An increased understanding of the phenomenon of 
polymorphism should enable pharmaceutical scien- 
tists to gain control over the crystallization process in 
order to selectively obtain the desired polymorph or 
suppress the growth of an undesired one. Phase 
changes during processing and scale-up are a prob- 
lem, which may be avoided by carefully designed 
initial small-scale studies. The availability of detailed 
structural data, combined with strategic design of 
substrates and additives, has led to significant ad- 
vances in the control over the polymorphs obtained 
in a particular crystallization [98]. With all the 
information available from these initial studies, it 
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should be possible to design and to select processing 
conditions which would give a desired polymorph 
and maintain the desired form throughout the various 
stages of drug processing and manufacture. 
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